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We present a review of the searches for very rare decays to muonic final states performed at 
^ LHCb using 1.0 fb _1 of pp collisions at 7 TeV centre of mass energy. Flavour changing neutral 

^ current processes, such as b9n — s- and B9 s —> are highly suppressed in the 

Oh Standard Model (SM). Such decays therefore allow contributions from new processes or new 

heavy particles to significantly modify the expected SM rates. Charged lepton flavour violating 
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processes, such as the neutrino-less T~ — > decay, have vanishingly small decay rates 

in the SM, but can be significantly enhanced in extended models. We report the latest results on 
these channels from LHCb. 
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1. Introduction 

The searches for B9v — > b9, — > }1 + }X~ ji + pi~ and %~ — > }i~ /i~ }X + performed by analysing 

1.0 fb _1 of pp collisions recorded in 201 1 with the LHCb detector [1] are presented. The three de- 
cay modes being purely muonic, the analyses have similarities. They all take advantage of the effi- 
cient muonic trigger lines [2], that for instance allow to trigger events with a muon of low transverse 
momentum (down to 0.5 GeV for di-muons), and reach an efficiency of 90% for di-muon signal. 
In addition, high rejections of combinatorial backgrounds are achieved using the excellent track- 
ing performances (e.g. the resolution on the momentum is ^ G [0.4,0.6] for p G [5; 100] GeV). 
Finally physical backgrounds, where hadrons are misidentified as muons, are reduced with the 
particles identification detectors that yield an efficiency for identifying muon of 97% for misiden- 
tification rate of pions (kaons) into muons below 3 (5)%. To avoid unconscious bias, the events for 
which the B9 ^ or % mass candidates he within three mass resolutions of the expected mean are kept 
hidden until all crucial choices for the analysis were made. 

2. Searches for B9^ ->• 

Within the Standard Model (SM), B9* — > are rare processes as they occur only via 

V s ) 

loop diagrams and are helicity suppressed. The amplitudes contributing to the branching frac- 
tion (38) can be expressed in terms of the scalar (cs), pseudoscalar (cp) and axial vector (ca) 
Wilson coefficients in a general approach [3]. Within the SM, c$ and cp contributions are neg- 
ligible while ca is calculated with a few percent accuracy [4] and leads [5] to a prediction of 
38{B Q ^ jU+aOsm = (3.2 ±0.2) x 10~ 9 and 38(B ^ a* + AOsm = (1.0 ±0.1) x lO" 10 . Mod- 
els beyond the SM could contribute to these Wilson coefficients and change significantly the 38. 
For instance, within the Minimal Supersymmetric SM (MSSM) in the large tan/3 approxima- 
tion [6], c^| SM o< tan 3 f5 /M\, where Ma denotes the pseudoscalar Higgs mass and tan j8 the ratio of 
Higgs vacuum expectation values. The LHCb experiment sets the most restrictive upper limits [7], 
38(B° S ^ At + AO< 1.4 x 10~ 8 and 38(B ^ jU + aO< 3-2 x 10~ 9 at 95% C.L. The present paper 
expounds the update of this analysis using the full 2011 data set. 1 

The strategy of the analysis is to derive the expected numbers of events of background and 
signal for a given 38 hypothesis and to compare these numbers to the observed ones with the CL S 
method [9]. The method provides CL S+ ^(CL^), a measure of the compatibility of the observed dis- 
tribution with the signal plus background (background only) hypothesis, and CL S = CL s+ ^/CL^ 
which is used to set upper limits on the 3$. 

For a given 38, the number of events is obtained by scaling the yields observed in three controls 
channels, B° — > J/V0, B°— > K + 7T~ and B+ — > J/i/^K + , by the ratio of the 38 hypothesis to the control 
channel 38. The scaling factor corrects also for the different efficiencies between signal and control 
channels and for different intial states using f s jf c \ measured at LHCb [10]. 

To improve the sensitivity of the analysis, a selection which reduces the backgrounds is ap- 
plied. Backgrounds stem from combinatorics where the two muons come from two different b- 
mesons and from peaking background from B9^ — > h + h'~ (h standing for K + or n + ) where hadrons 

'Since then, the analysis has been updated leading to the first evidence for the B^— !• /l + /i~decay [8], 
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Figure 1: Distribution of selected candidates (black points) in the (left) — >• jU + jU~ and (right) B°— > jU + jU~ 

mass window for BDT>0.5, and expectations for, from the top, B|? — > SM signal (gray), combinatorial 

background (light gray), B9 , — > h + h'~ background (black), and cross-feed of the two modes (dark gray). The 

hatched area depicts the uncertainty on the sum of the expected contributions. 



are identified as muons. The first type of background is reduced by cutting on topological and kine- 
matical variables and on a combination of them obtained with a Boosted Decision Tree (BDT). This 
selection has an efficiency similar on the signal, normalisation and control channels. The peaking 
background is reduced using information from the particle identification detectors. 

After this selection, events are classified in bins of the di-muon invariant mass and bins of 
a topological variable built with a second BDT. In each bin the expected numbers of signal and 
background events are derived. The signal mass shape is assumed to be a Crystal Ball shape, 
where the mean and the resolution are obtained on data while the transition point is derived from 
simulations. The BDT shape is obtained on data by extracting with a fit to the mass distribution, 
in each BDT bin, the yields of B^ — > h + h'~. For the combinatorial background, the mass and the 
BDT shape are obtained simultaneously by interpolating, in each BDT bin, the mass side-bands 
into the signal regions with an exponential function. Finally, the peaking background BDT shape 
is assumed to be the same as the signal one (correcting for trigger bias) and the mass shape is taken 
from simulations. The total number of peaking backgrounds is derived by a data driven method. 
The observed data and the expected distributions for the SM signal and background components 
are shown in Figure 1 . 

The comparison of the distributions of observed events and expected background events results 
in a p-value (1 - CL b ) of 18 % (60 %) for the B° -> /i+jU" (B° -» jU+ju") decay, where the CL b 
values are those corresponding to CL s+b = 0.5. No excess being observed, uppers limits are set and 
reported in Table 1. The expected and observed CL S values are shown in Fig. 2 for the B ( s '— > [l + pL~ 
and B° — > pi + pL~ channels, each as a function of the assumed 2$. 

Assuming the same sensitivity for the analysis, prospect for the discovery of the B° — > }JL + pL~ 
decay can be derived. Figure 3 shows the value of ^(B® ~ > M + M ) f° r which a 3 a evidence of 
the decay could be obtained as a function of the integrated luminosity. This extrapolation takes into 
account the increase of energy at the center of mass to 8 TeV and the already observed data pattern 
for 2011. 
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Figure 2: CL S as a function of the assumed SS for (left) — > n + /j.~and (right) B°— > j(i + ^~decays. The 
long dashed black curves are the medians of the expected CLs distributions for B|? — > if background 

and SM signal were observed, and for B° — > n + jj,~, if background only was observed. The yellow areas 
cover, for each 34% of the expected CLs distribution on each side of its median. The solid blue curves 
are the observed CLs. The upper limits at 90% (95 %) C.L. are indicated by the dotted (solid) horizontal 
lines in red (dark gray) for the observation and in gray for the expectation. 



Mode 
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B°^ju 
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at 90%C.L. 


at 95 % C.L. 


at 90 % C.L. 


at 95 % C.L. 


Exp. bkg+SM 


6.3 x 1(T 9 


7.2 x 1(T 9 






Exp. bkg 


2.8 x 1(T 9 


3.4 x 1(T 9 


0.91 x 10~ 9 


1.1 x 10~ 9 


Observed 


3.8 x 1(T 9 


4.5 x 1(T 9 


0.81 x 10~ 9 


1.0 x 10~ 9 



Table 1: Expected and observed limits on ^(B? > -> n + n ). 
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Figure 3: 88 at which the B^ — > ju + jtt~ decay would be observed at 3d as a function of the integrated 
luminosity. The black plain line is the median of the 38 distribution giving a 3(7 observation. The yellow 
areas cover, for each integrated luminosity, 34% of the SB distribution on each side of its median. The red 
horizontal band represents the time integrated SM prediction for ^(B^ — > [11]. 

3. Searches for B9 , ->■ /i + fi~ /i + pT 

The decays of B9 ^ mesons to four muons are strongly suppressed in the SM and, for B°, 
mainly occur via the resonant mode B ( s ' — y J/i//(/i + /x _ )<I>( i u + /x _ ) for which the 88 is estimated 
to be (2.3 ±0.9) x 10 [12]. In the SM, the $M of the non-resonant modes are predicted not 
to exceed 10~ 10 [13], but could receive non negligible contributions from models beyond the SM. 
These decay modes have therefore potential for discovering physics beyond the SM and reassessing 
the hyper-CP anomaly [14]. 
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Figure 4: Observed invariant mass distribution of the B9, 
area is the B°(B°) mass window. 



• pL + pL pL + pi candidates. The blue (red) dashed 



The searches are implemented with a cut and count analysis. Similarly to the B" s j — > }l + }X 

analysis, a control channel, B°— > J/vA(/x/i)K*°(K + 7r _ ), is used to derived the number of signal 
events corresponding to a S§ hypothesis. The selection is designed to be similar for the signal and 
the normalisation channel. The four muon candidates are required to be made of tracks with a 
good particle identification and to have an impact parameter significance (IPS) with respect to the 
primary vertex larger than 16. The candidate should also form a vertex with a % 2 per degree of 
freedom smaller than 6 and have an IPS smaller than 9. Resonant candidates are vetoed when their 
di-muons invariant masses lie around the masses of the J/y and <1>. 

The only relevant source of background is the combinatorial one. The expected number of 
these events is obtained by interpolating the side -bands and yields to 0.38+QJ7 (0.30^q event s 
in the B° (B°) mass window. The observed pattern is shown in Figure 4 and is compatible with the 
background expectation. Using the CLs method [9] the first limits on these two processes are set 
to^(B°->- H+/X-H+IX-) < 1.3 x 1(T 8 and^(B°^jU + j"~Ai + AO <5.4x 1(T 9 at 95% C.L. . 



4. Search for x ->■ \i ji \i + 

In the SM the lepton flavour violation (LFV) decay %~ — > pL~ fi~ jJ. + is almost completely 
suppressed, while many theories beyond the SM [15] predict an enhanced LFV and in particular in 
the T sector, raising in some cases g${x~ — > to an experimentally observable level. 

The analysis strategy is similar to the one used for the searches of B9 ^ — > jU + jU '. The main 
background processes consist of long decay chains of heavy mesons with three real or misidentified 
muons in the final state. The dominant contribution in these processes is coming from D.T — > 
tj(/i + \i~ 'y)n~v : fi, and is treated separately. 

After a loose selection based on the kinematics of the decay, the three muon candidates 
are classified in bins of invariant mass and of two multivariate operators combining information 
of the topology of the candidates and of the particle identification of the tracks used to form 
them. The three variables line shapes are determined for the signal with data driven methods 
and for the background by interpolating from the mass side-bands or with simulated data for 
DjT — > r\(}x + pL~ y)jx~v ^. The signal yield corresponding to a SB hypothesis is obtained by nor- 
malising to D.T — > <£>(h + ij.~)k~. 
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Figure 5: (left) Observed invariant mass distribution of the T~ — > pL~ pL~ ji + candidates in the most signal 
like merged bins. The solid line is the fit to the data of the sum of the combinatorial exponential contribution 
(dashed red line) and of the D~ — > r\(fi + fi~ y)fi~v ^ one (dashed green line), (right) CL S as a function 
of the assumed SS. The long dashed blue curves are the medians of the expected CLs distributions for 
T~ — > fl~ fl~ }l + if background only was observed. The yellow (green) areas cover, for each SS, 34% (48%) 
of the expected CLs distribution on each side of its median. The solid black curves are the observed CLs. 

The observed data together with the background expectations are shown in Figure 5. The com- 
patibility of the data to the background or background plus signal hypothesis are computed with the 
CL S method [9]. The observed data are compatible with the background only hypothesis. Figure 5 
shows the CL S as a function of the assumed An upper limit is set at ^(t~ — > pL~ pL~ pL + ) < 
7.8(6.3) x 10~ 8 at 95 (90)% C.L. which is the first limit of a LFV X decay set at a hadron collider. 
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